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Abstract 

In this paper, the deployment of an unmanned aerial vehicle (UAV) as a flying base station used to 
provide on the fly wireless communications to a given geographical area is analyzed. In particular, the 
co-existence between the UAV, that is transmitting data in the downlink, and an underlaid device-to- 
device (D2D) communication network is considered. For this model, a tractable analytical framework 
for the coverage and rate analysis is derived. Two scenarios are considered: a static UAV and a mobile 
UAV. In the first scenario, the average coverage probability and the system sum-rate for the users in the 
area are derived as a function of the UAV altitude and the number of D2D users. In the second scenario, 
using the disk covering problem, the minimum number of stop points that the UAV needs to visit in 
order to completely cover the area is computed. Furthermore, considering multiple retransmissions for 
the UAV and D2D users, the overall outage probability of the D2D users is derived. Simulation and 
analytical results show that, depending on the density of D2D users, optimal values for the UAV altitude 
exist for which the system sum-rate and the coverage probability are maximized. Moreover, our results 
also show that, by enabling the UAV to intelligently move over the target area, the total required transmit 
power of UAV while covering the entire area, is minimized. Finally, in order to provide a full coverage 
for the area of interest, the tradeoff between the coverage and delay, in terms of the number of stop 
points, is discussed. 


I. Introduction 

The use of unmanned aerial vehieles (UAVs) as flying base stations that ean boost the eapaeity 
and coverage of existing wireless networks has recently attracted significant attention [0 and 
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Q. One key feature of a UAV that ean potentially lead to the eoverage and rate enhaneement is 
having line-of-sight (LoS) eonneetions towards the users. Moreover, owing to their agility and 
mobility, UAVs ean be quiekly and effieiently deployed to support eellular networks and enhanee 
their quality-of-serviee (QoS). On the one hand, UAV-based aerial base stations ean be deployed 
to enhanee the wireless eapaeity and eoverage at temporary events or hotspots sueh as sport 
stadiums and outdoor events. On the other hand, they ean be used in publie safety seenarios 
to support disaster relief aetivities and to enable eommunieations when eonventional terrestrial 
networks are damaged Q . Another important applieation of UAVs is in the Internet of Things 
(loT) in whieh the deviees have small transmit power and may not be able to eommunieate over 
a long range. In this ease, a UAV ean provide a means to eolleet the loT data from one deviee 
and transmit it to the intended reeeiver 0 and 0 . Last but not least, in regions or eountries in 
whieh building a eomplete eellular infrastrueture is very expensive, deploying UAVs is highly 
benefieial as it removes the need for towers and eables. In order to reap the benefits of UAV 
deployments for eommunieation purposes, one must address a number of teehnieal ehallenges 
that inelude performanee analysis, ehannel modeling, optimal deployment, resouree management, 
and energy effieieney, among others [|5|-p^. 

The most signifieant existing body of work on UAV eommunieations foeuses on air-to-ground 
ehannel modeling 0-0. For instanee, in 0 and 0, the probability of line of sight (LoS) for 
air-to-ground eommunieation as a funetion of the elevation angle and average height of buildings 
in a dense urban area was derived. The air-to-ground path loss model has been further studied 
in 0 and 0. As diseussed in 0, due to path loss and shadowing, the eharaeteristies of the 
air-to-ground ehannel are shown to depend on the height of the aerial base stations. 

To address the UAV deployment ehallenge, the authors in 0 derived the optimal altitude 
enabling a single, statie UAV to aehieve a maximum eoverage radius. However, in this work, 
the authors simply defined a deterministie eoverage by eomparing the path loss with a speeified 
threshold and did not eonsider the eoverage probability. The work in p0| extends the results 
of 0 to the ease of two UAVs while eonsidering interferenee between the UAVs. In pT| , 
the authors studied the optimal plaeement of UAVs for publie safety eommunieations in order 
to enhanee the eoverage performanee. However, the results presented in GD are based on 
simulations and there is no signifieant analytieal analysis. Moreover, the use of UAVs for 


supplementing existing eellular infrastrueture was diseussed in [12| whieh provides a general 
view of praetieal eonsiderations for integrating UAVs with eellular networks. The work in p0 
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considered the use of UAVs to eompensate for the eell overload and outage in eellular networks. 


However, [ 121 does not provide any analysis on the eoverage performanee of UAVs and their 
optimal deployment methods. In p^, the authors investigated how to optimally move UAVs for 


improving eonneetivity of ad hoe networks. However, [14| only foeused on an ad-hoe network 


and assumed that the UAV have eomplete information about the loeation of nodes. In [I5|, 
eonsidering statie ground users, the optimal trajeetory and heading of UAVs equipped with 
multiple antennas for ground to air uplink seenario was derived. The work in proposed a 
power effieient deployment and eell assoeiation for multiple UAVs in downlink transmissions. 

For seenarios in whieh there is limited or no infrastrueture support, beyond the use of 
UAVs, there has been eonsiderable reeent works that study the use of direet deviee-to-deviee 
(D2D) eommunieations between wireless users over the lieensed speetrum (n). Sueh D2D 
eommunieations has been shown improve eoverage and eapaeity of existing wireless networks, 
sueh as eellular systems. In partieular, in hotspot areas or publie safety seenarios, D2D will 
allow users to eommunieate direetly with one another without signifieant infrastrueture. D2D 
eommunieations are typieally deployed using underlaid transmission links whieh reuse existing 


lieensed speetrum resourees [18|. Therefore, deploying a UAV over a speetrum band that must 
be shared with an underlaid D2D network will introduee important interferenee management 
ehallenges. In the literature, there are some studies on the eoexistenee of the underlaid D2D and 


eellular eommunieations with a single base station [ 191. Furthermore, the authors in [20| exploited 
the interplay between the massive MEMO and underlaid D2D eommunieations for a single eell 


ease. The authors in [211 extended the previous work on the D2D/massive MIMO eoexistenee to 
the multi-eell seenario. However, none of theses prior works studied the eoexistenee of UAVs and 
underlaid D2D eommunieations. In partieular, a eomprehensive analytieal analysis to evaluate 
this eoexistenee in terms of different performanee metries, sueh as eoverage and rate, is laeking 
in the eurrent state-of-the-art @, pAl , [[T^-[|2T|. 

Compared to the previous studies on the eoexistenee of D2D and eellular networks sueh as 


and [CT|, the presenee of an aerial UAV base station along with D2D links introduees new 
ehallenges. First, the ehannel modeling between the UAV and ground users will no longer be 
a elassieal fading ehannel, instead, it will be based on probabilistie LoS and NLoS links [[^, 
0, while the ehannel between a base station and the users will still follow a Rayleigh fading 
model. Seeond, unlike eonventional, fixed base stations, the height of a UAVs is adjustable 
and this impaets the ehannel eharaeteristies and the eoverage performanee. Third, the potential 
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mobility of a UAV introduces new dimensions to the problem and the impaet of sueh mobility 
on D2D and network performanee must be analyzed. The prior studies on UAVs sueh as 0-G3 
have not addressed the third ehallenge. More speeifieally, the interplay between UAVs and D2D 
eommunieations and the existing ehallenges and tradeoffs have not been investigated in these 
literature. To our best knowledge, this paper will provide the first comprehensive fundamental 
analysis on the performanee of UAV eommunieation in the presenee of underlaid D2D links. 

The main eontribution of this paper is to analyze the eoverage and rate performanee of 
UAV-based wireless eommunieation in the presenee of underlaid D2D eommunieation links. 
In partieular, we eonsider a network in whieh a single UAV must provide downlink transmission 
support to a number of users within a given area. In this area, a subset of the deviees is also 
engaged in D2D transmissions that operate in an underlay fashion over the UAV’s transmission. 
We eonsider two types of users, namely downlink users (DUs) whieh reeeive data from the 
UAV, and D2D users whieh eommunieate direetly with one another. Here, the UAV must 
eommunieate with the DUs while taking into aeeount the potential interferenee stemming from 
the underlaid D2D transmissions. For this network, we analyze two key eases: statie UAV and 
mobile UAV. Using tools from stoehastie geometry, for both seenarios, we derive the average 
downlink eoverage probabilities for DUs and D2D users and we analyze the impaet of the UAV 
altitude and density of the D2D users on the overall performanee. For the statie ease, we find the 
optimal values for the UAV altitude whieh leads to a maximum eoverage probability for DUs. 
In addition, eonsidering both DUs and D2D users, an optimal altitude whieh maximizes the 
system sum-rate is eomputed. Our results demonstrate that the optimal UAV altitude deereases 
as the density of D2D users inereases. The results show that a maximum system sum-rate ean 
be aehieved if the UAV altitude is appropriately adjusted based on the D2D users’ density. 
Furthermore, for a given UAV altitude, we show that an optimal value for the number of D2D 
users that maximizes the system sum-rate exists. 

For the mobile UAV ease, we assume that the UAV ean travel over the area while stopping 
at some given loeations in order to serve the downlink users. Using the disk eovering problem, 
we find a minimum number of stop points that the UAV needs to to eompletely eover the area. 
This ean be interpreted as the fastest way to eover the whole area with a minimum required 
transmit power. In addition, we analyze the tradeoff between the number of stop points, whieh 
is eonsidered as delay here, and the eoverage probability for the downlink users. Moreover, 
eonsidering retransmissions at different time instanees, we derive the overall outage probability 
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for the D2D communications. We show that, in order to enhance the coverage for DUs, the UAV 
should stop in more locations over the target area which can, in turn, lead an increased delay for 
DUs and higher outage probability for D2D users. For example, our results show that for a given 
density of D2D users, to increase the DU coverage probability from 0.4 to 0.7, the number of 
stop points should be increased from 5 to 23. Furthermore, the number of stop points is shown 
to significantly depend on the number of D2D users. For instance, if the average number of D2D 
users in the area increases from 50 to 100, in order to maintain the DUs’ coverage requirement, 
the number of stop points should be increased from 20 to 55. 

The rest of this paper is organized as follows. Section II presents the system model and 
describes the air-to-ground channel model. In Section III, coverage probabilities for DUs and 
D2D users are provided for a single static UAV. Section IV presents the performance evaluation 
for one mobile UAV which is used to provide full coverage for the target area. Section V presents 
the simulation results while Section VI draws some conclusions. 


II. System Model 

Consider a circular area with a radius Rc in which a number of wireless users are deployed. 
In this area, a UAV (at low altitude platform) is deployed to act as a flying base station and 
serve a subset of those users. In this network, the users are divided into two groups: downlink 
users located uniformly in the cell with density Xdu (number of users per m^), and D2D users 
whose distribution follows homogeneous Poisson point processes (PPP) $b with density 
of Xd (number of D2D pairs per m^). Note that, the average number of users in a given area 
is equal to the density of the users multiplied by the size of the area. Here, we focus on 
the downlink scenario for the UAV and we assume that the D2D users communicate in an 
underlay fashion. Furthermore, we assume that a D2D receiver connects to its corresponding 


D2D transmitter pair located at a fixed distance away from it in an isotropic direction p9|. 
Therefore, the received signals at the D2D receiver include the desired signal from the D2D 
transmitter pair and interference from the UAV and other D2D transmitters. A downlink user, 
on the other hand, receives the desired signal from the UAV but it also experiences interference 
from all the D2D transmitters. It should be noted that, in our model, the UAV provides service 
for downlink users (DUs) located inside a given, finite area with radius Rc- Nonetheless, we 
assume that the D2D users are spatiality distributed according to a PPP over an infinite area. 
In other words, each user receives interference from an infinite number of D2D transmitters. 
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Figure 1: Network model ineluding a UAV, downlink users and D2D. 


This is a typieal assumption in PPP analysis whieh ensures that, the average amount of reeeived 


interferenee from D2D transmitters does not depend on the loeation of the users [20|, [231, and 


[24|. 


The signal to interferenee plus noise ratio (SINK) expression for a D2D reeeiver is: 

Pr,d 


Id = 


P + Iu + N^ 


( 1 ) 


where Pr^ is the reeeived signal power from the D2D transmitter, is the total interferenee 
from other D2D users, is the interferenee from the UAV, and N is the noise power. Moreover, 
we have: 


Pr,d = Pddo “‘"S'O, 

(2) 


(3) 



Id = J2Pddr"'^g^, 

(4) 


i 


where the index z = 0 is used for the seleeted D2D transmitter/reeeiver pair, qq and Qi are, 
respeetively, the ehannel gains between a D2D reeeiver and its eorresponding D2D transmitter, 
and the interfering D2D transmitters. For the D2D transmission, we assume a Rayleigh fading 
ehannel model [ [T^ , [ |2T| and [25|. Pd is the D2D transmit power whieh is assumed to be fixed 
and equal for all the users, di is the distanee between a D2D reeeiver and the z^^ D2D transmitter, 
do is the fixed distanee between the D2D reeeiver and transmitter of the seleeted D2D pair, and 
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ad is the path loss exponent between D2D users. Note that the reeeived signal powers as well 
as the noise power are normalized by a path loss eoefficient. 

The SINK expression for a DU user that conneets to the UAV is: 


lu 


Id + N^ 


(5) 


where is the received signal power from the UAV and Id is the total interference power 
from D2D transmitters. Finally, the SINR-based coverage probability for the downlink users and 
the D2D users is given by: 


Pcov,du{/3) = P [7« > , 

(6) 

Pco.,dW) = P [7a > /S], 

(7) 


where 7 „ and 7 ^ are, respectively, the SINK values at the location of the downlink users and 
the D2D users, and (3 is the SINK threshold. 


A. Air-to-ground channel model 

As discussed in [|^ and the ground receiver receives three groups of signals including 
LoS, strong reflected non-line-of-sight (NLoS) signals, and multiple reflected components which 
cause multipath fading. These groups can be considered separately with different probabilities 
of occurrence as shown in Q and Q. Typically, it is assumed that the received signal is 
categorized in only one of those groups [[^. Each group has a specific probability of occurrence 
which is a function of environment, density and height of buildings, and elevation angle. Note 
that the probability of having the multipath fading is significantly lower than the LoS and NLoS 
groups 0 . Therefore, the impact of small scale fading can be neglected in this case Q. One 
common approach to modeling air-to-ground propagation channel is to consider LoS and NLoS 
components along with their occurrence probabilities separately as shown in Q and Q. Note that 
for NLoS connections due to the shadowing effect and the reflection of signals from obstacles, 
path loss is higher than in LoS. Hence, in addition to the free space propagation loss, different 
excessive path loss values are assigned to LoS and NLoS links. Depending on the LoS or NLoS 
connection between the user and UAV, the received signal power at the user location is given 
by 0 : 



— 0!.u 


LoS connection, 
NLoS connection, 


( 8 ) 


P 



where is the UAV transmit power, |X„| is the distance between a generic user and the UAV, 
au is the path loss exponent over the user-UAV link, and r] is an additional attenuation factor due 
to the NLoS connection. Here, the probability of LoS connection depends on the environment, 
density and height of buildings, the location of the user and the UAV, and the elevation angle 
between the user and the UAV. The LoS probability can be expressed as follows 0: 


1 


(9) 


l + Uexp(-5[0-C'])’ 


where C and B are constant values which depend on the environment (rural, urban, dense urban. 



or others) and 6 is the elevation angle. Clearly, 9 = ^ x sin“^ ~ \/JP + and 

also, probability of NLoS is PnloS = 1 ~ PloS- 

As observed from (9), the LoS probability increases as the elevation angle between the user and 
UAV increases. 

Given this model, we will consider two scenarios: a static UAV and a mobile UAV. For each 
scenario, we will derive the coverage probabilities and average rate for DUs and D2D users. 
Once those metrics are derived, considering the D2D users density, we obtain optimal values 
for the UAV altitude that maximize the coverage probability and average rate. 


III. Network with a Static UAV 


In this section, we evaluate the coverage performance of the scenario in which one UAV 
located at the altitude of h in the center of the area to serve the downlink users in the presence 
of underlaid D2D communications. It can be shown that, for a uniform distribution of users over 
the given area, palcing the UAV in the center of the area can maximize the coverage probability 
of the downlink users. 

A. Coverage probability for D2D users 

Consider a D2D receiver located at (r, p), where r and p are the radius and angle in a polar 
coordinate system assuming that the UAV is located at the center of the area of interest. The 
distance between the D2D transmitter and its corresponding receiver is fixed and it is denoted 
by do- In this case, for underlaid D2D communication, the coverage probability for the D2D 
users can be derived as follows: 
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Theorem 1. The coverage probability for a D2D receiver, at the location (r, 0), connecting to 
its D2D transmitter located at a distance do away from it, is given by: 


-Pcov,d(r, = exp 


—27r' 




X I^PLosexp 

where |X„| = + r^. 


adSm{2Tr/ad) Pd J 

-Pd^-^PfX, 


Pd 


+ ^NLoS exp 


Pd 


, (10) 


Proof See Appendix A. ■ 

From this theorem, we ean make several key observations. First, eonsidering the faet that the 
UAV ereates interferenee on the D2D users, inereasing the UAV altitude to inerease its distanee 
from the D2D users does not neeessarily reduees the interferenee on the D2D users. As will 
be shown later by numerieal simulations, by inereasing the UAV altitude the D2D eoverage 
probability deereases first, and then inereases. This is due to the faet that, eonsidering (9) and 
(10), although inereasing the UAV altitude inereases the path loss term, it also leads to a LoS 
probability. In general, the D2D users prefer to have the NLoS view towards the UAV and have 
a maximum distanee from it, however, these two objeetives eonfliets with eaeh other. Seeond, 
inereasing the D2D transmit power {Pd), always enhanees the D2D eoverage probability, even 
in an interferenee limited seenario where noise is ignored. Typieally, in the interferenee limited 
seenarios, inereasing the transmit power of the D2D users does not improve the eoverage per- 
formanee due to the inereased interferenee from other D2D transmitters. Aeeording to Theorem 
1, although in the interferenee limited seenario {N = 0) the first multiplying term in (10) is 
independent of Pd due to the interferenee from D2D transmitters, the seeond term is an inereasing 
funetion of Pd- Finally, the D2D eoverage probability in (10) deereases when the UAV transmit 
power inereases. To eope with this situation, the D2D users ean inerease their transmit power or 
reduee the fixed distanee parameter (do)- In addition, deereasing the D2D user density improves 
the eoverage probability due to deereasing the interferenee. Note that the result presented in 
Theorem 1 eorresponds to the eoverage probability for a D2D user loeated at (r, p). To eompute 
the average eoverage probability in the eell, we eonsider a uniform distribution of users over 
the area with f{r,(p) = 0<r<Pc,0<(p< 27ij^ where Rc is the radius of the 

'Note that the number of users has a Poisson distribution but their location follows the uniform distribution over the area 
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desired eireular area. Then, we eompute the average over the desired area. The average eoverage 
probability for D2D users will be: 

PcovAl^) = /3)] 


Rc 


= exp 


= exp 


-27T^X,(3P^‘^dl /3CiV \ ^ 

arfsin(27r/ad) Pd ) J 


-27i^Xdf3P^‘^dl _ ^d^\ 


Rc 


arfsin(27r/ad 


Pd J 


E 


lu 


exp(- - -) 

^d 

-/5C4; 


exp( 


Pd 


/(r, ip)drdip 


2r , 


(11) 


From (11), we ean see that the average eoverage probability for D2D users inereases as the 
size of the area, Rc, inereases. In faet, when the UAV serves a larger area, the average distanee 
of D2D users from the UAV inereases and on the average they reeeive lower interferenee from 
it. Next, we provide a speeial ease for (11) in whieh the UAV has a very high altitude or very 
small transmit power. 

Remark 1. For = d or h ^ oo, the average coverage probability for the D2D users is 

'-2TT^Xd0^‘^dl 


simplified to [24]: 


PcovA/^) = exp 


adSin(27r/Q;d) 


Pd 


( 12 ) 


Note that, (12) corresponds to the coverage probability in overlay D2D communication in which 
there is no interference between the UAV and the D2D transmitters. It should be noted that, this 


result is also related to the success probability in a bipolar ad hoc network [24], 


B. Coverage Probability for Downlink Users 

Here, we first derive the bound and lower bound for the downlink users’ eoverage probability. 

Theorem 2. The lower bound and upper bound of the average coverage probability for DUs in 
the area of interest is given by: 


Rc 


PL,AAh)= / PhoB{r,h)Li 


0 

Rc 


Pu\X,, 


+ / -PnLos(^) 


pPuAc, 


/S 


2r 

-ivl dr 


_ iV ) —dr 


( 13 ) 


Rc 


^cov,du(/5,^)= PLos{r,h)Ui 


Pu\X,, 




N 1 |dr 


0 




















( 14 ) 


+ J PNLosir,h)Ui 

0 

where f3N < and for any T > 0, 

\ 2 y Pd J 


'nPu\Xn 


/S 


T 


N ) —dr 
I ar, 


-"^lOLd 


exp I -7rXd[ — ) r(l + 2/ad) ) , (15) 


Ui{T) = exp I -vrAdl — 


— ‘i/Old 


r(l + 2 / ad) I . 


(16) 


Also, r(t) = J ^dx is the gamma function |26|. 

0 

Proof. See Appendix B. 


From Theorem 2, we can first see that, for T » Pd, given that e ^ 1 — x when x —0, 

we have Ui{T) = Lj{T) 1 — r(l + "^joid). This means that the lower bound 

and upper bound become tighter for lower transmit power of D2D users. Moreover, from (15) 
and (16), when \d —)■ cxd, the number of D2D users tends to infinity and Ui = Lj = 0. 
Consequently, the downlink users experience an infinite interference from the D2D users which 
results in Pcov.du = 0. 

Furthermore, considering (9), (13), and (14), we can see that increasing the UAV altitude 
Qi), can enhance the LoS probability and the coverage probability. On the other hand, due to 
increasing \Xu\, Lj and Uj decrease, and hence the coverage probability for downlink users 
decreases. Therefore, in order to achieve the maximum coverage, the altitude of the UAV should 
be carefully adjusted. 

As per Theorem 2, increasing Rc decreases the average coverage probability for the downlink 
users. However, higher Rc results in a higher D2D average coverage probability. Moreover, 
the average coverage probability for downlink users decreases as the density of the D2D users 
increases. In this case, to improve the DUs coverage performance, one must increase or 
reduce Rc. Next, we derive the DU coverage probability in the absence of the D2D users. 


Proposition 1. For low density and transmit power of D2D users, the interference from D2D 
users is negligible compared to the UAV, then, the exact average coverage probability for the 
downlink users can be expressed by: 


[P) 


■ \i Pf, 1 

Ad 
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Proof. For a DU located at (r, (^), the eoverage probability in absenee of D2D users beeomes 

Pcov,du(u /3) = p [ 7 . >13] = PLos(r)P [7u > /3|LoS] + PNLos(r)P [ 7 . > /3|NLoS] 


= Plos(^^)1 


r < 


P, 




l/o„ 


/3ivy 


+ PnLos(’" )1 


■<'s) 


l/o„ 


( 18 ) 


The average eoverage probability is eomputed by taking the average of Pcov,du{r, (f, f3) over the 
eell with the radius Pc- 

PcoN,du{j'i Pi (3') [PcoVjduiXi Pi / 3 )] 




2r 2r 

PLos{r)—dr + J PNLos{r)—dr. (19) 


Proposition 1 gives the exaet expression for the downlink users’ eoverage probability when 
the interferenee from D2D users, due to their low density and low transmit power, is negligible 
eompared to the UAV. Therefore, the result in Proposition 1 shows the maximum aehievable 
eoverage performanee for downlink users when the reeeived signal from the UAV is dominant 
eompared to the interferenee from the D2D transmitters. 


C. System sum-rate 

Now, we investigate the average aehievable rates for the DUs and D2D users whieh ean be 


expressed as in p0| : 

Cdu = W\og2 (1 + /?) Pcov,dn (/3)i (20) 

U, = lUlog 2 (l + /?)Pcov,d(/9), (21) 

where W is the transmission bandwidth. Considering the whole DUs and D2D users in the eell, 
the system sum-rate, Usum, can be derived as a funetion of the eoverage probabilities and the 


number of users as follows: 


C’sum Pc ^3\duCdu P Pc ^3\(iCd- 


( 22 ) 


Assuming /i = ^, we have 

■^d 

Gsum = ArfPcV [pPcov,du{l3) + Pcov,<i(/5)] fFlog 2 (l + /3), (23) 

where Pc^vrA^ and RfirSdu are the number of DUs and D2D users in the target area respeetively. 
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From (30), observe that, on the one hand, Csum is direetly proportional to A^, but on the other 
hand, it depends on the eoverage probabilities of DUs and D2D users whieh both are deereasing 
funetions of D2D user density. Therefore, in general inereasing does not neeessarily enhanee 
the rate. Note that, eonsidering (11), (13), (14), and (23), for both A^ —)■ 0 and Xd ^ oo eases 
the system sum-rate tend to zero. Henee, there is an optimum value for A^ that maximizes Csum- 
Aeeording to (23), Csum is a funetion of the eoverage probability and a logarithmie funetion 
of the threshold ((3). The former is a deereasing funetion of (3 whereas the latter is an inereasing 
funetion of (3. In other words, although inereasing the threshold is desirable for the rate due 
to inereasing the logarithmie funetion, it also reduees the eoverage probability. Therefore, in 
order to aehieve a maximum rate, a proper value for the threshold ean be adopted. It should be 
noted that, the SINK threshold, j3, is typieally fixed and eannot be set lower than the reeeiver 
sensitivity. However, the analysis of different values of f3 brings value in order to understand how 
one eould ehange the SINK threshold value (in the future) through proper resouree alloeation or 
just system design (ehange in the number of users, ete). 


IV. Network with a Mobile UAV 

Now, we assume that the UAV ean move around the area of radius in order to provide 
eoverage for all the downlink users in the target area. In partieular, we eonsider a UAV that 
moves over the target area and only transmits at a given geographieal loeation (area) whieh we 
hereinafter refer to as “stop points”. Eaeh stop point represents a loeation over whieh the UAV 
stops and serves the present downlink users. Here, our first goal is to minimize the number of 
stop points (denoted by M) and determine their optimal loeation. Note that, as the UAV moves, 
it ean have a different ehannel to a user at different time instanees. The objeetive of the UAV is 
to eover the entire area and ensure that the eoverage requirements for all DUs are satisfied with 
a minimum UAV transmit power and minimum number of stop points. In other words, we find 
the minimum number and loeation stop points for the UAV to eompletely eover the area. We 


model this problem by exploiting the so-ealled disk covering problem [27|. In the disk eovering 
problem, given a unit disk, the objeetive is to find the smallest radius required for M equal 
smaller disks to eompletely eover the unit disk. In the dual form of the problem, for a given 
radius of small disks, the minimum number of disks required to eover the unit disk is found. 

In Figure 2, we provide an illustrative example to show the mapping between the mobile 
UAV eommunieation problem and the disk eovering problem. In this figure, the eenter of small 
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disks can be considered as the location of stop points and the radius of the disk is the coverage 
radius of the UAV. Using the disk covering problem analysis, in Table I, we present, for different 
number of stop points, the minimum required coverage radius of a UAV for completely covering 


the target area [27|, [28|. Thereby, using the dual disk covering problem, for a given maximum 
coverage radius of a UAV, we can find the minimum number of stop points for covering the 
entire area. The detailed steps for finding the minimum number of stop points are provided next. 

First, we compute the coverage radius of the UAV based on the minimum requirement for 
the DU coverage probability. The coverage radius is defined as the maximum radius within 
which the coverage probability for all DUs (located inside the coverage range) is greater than 
a specified threshold, e. In this case, the UAV satisfies the coverage requirement of each DU 
which is inside its coverage range. The maximum coverage radius for the UAV at an altitude h 
transmitting with a power will be given by: 

Rm = max{i?|Pcov,dn(/3, R) > Pu, h} = P~l e), (24) 


where e is the threshold for the average coverage probability in the cell (area covered by the 
UAV). Note that, a user is considered to be in coverage if it is in the coverage range of the UAV. 
The minimum required number of stop points for the full coverage is: 


L = 

Pcov,duij') ^^ 

where M represents the number of stop points, the second condition guarantees that the area is 
completely covered by the UAV, and L is the minimum value for the number of stop points if 
the following condition holds: 

Pmin,L < Pm < Pmin,L-l min{M} = L. (26) 

By using Table I, we see that, Pmin,L-i and Pmin,L are, respectively, the minimum radius required 
to cover the entire target area with L — 1 and L disks. After finding the minimum M, we can 
reduce the UAV transmission power such that the coverage radius decreases to the minimum 
required radius (Pmin,^)- In this way, the UAV transmit power is minimized. Thus we have 

= argmin{P“^ ^^(/3, e) = R^,r,,L\h}, (27) 

Pu 

where Pu,min is the minimum UAV transmit power. Thereby, the minimum number of stop points 
leads to a full coverage at a minimum time with a minimum required transmit power. 

In summary, the proposed UAV deployment method that leads to the complete coverage with 
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Target area 



Figure 2: Five disks eovering problem. 


Table I: Number and radii of disks in the eovering problem. 


Number of stop points 

Minimum required coverage radius 

M = 1,2 

Rc 

CO 

9 

II 

9 -^C 

to 

II 

0.61i?e 

M = 6 

0.556i?c 

M = 7 

0.5i?e 

M = 8 

0.437i?c 

M = 9 

0.422i?c 

M = 10 

0.398i?, 

M = 11 

0.38i?e 

M = 12 

0.361i?c 


a minimum time and transmission power proeeeds as follows. First, depending on the parameters 
of the problem sueh as density of users and threshold, we eompute the maximum eoverage radius 
of a UAV at the optimal altitude that ean serve the DUs. Seeond, eonsidering the size of target 
area, using the disk eovering problem, we find the minimum required number of transmission 
points along with the eoverage radius at eaeh point. Third, we reduee the transmission power 
of UAV sueh that its maximum eoverage radius beeomes equal to the required eoverage radius 
found in the previous step. Using the proposed method, the target area ean be eompletely eovered 
by the UAV with a minimum required transmit power and minimum number of stop points. Next, 
we investigate the impaet of the number of stop points on the full eoverage time of the downlink 
users, and the overall outage probability of the D2D users. 

We eonsider the network during M time instanees in whieh the UAV and D2D users will 
exeeute M retransmissions. Note that, our system model eonsiders the downlink, therefore, 
the retransmissions are essentially from the UAV to the DUs, and from D2D transmitters to 
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corresponding receivers. The moving UAV satisfies the coverage requirements of the downlink 
users in M retransmissions from different locations. Clearly, as the number of stop points (M) 
increases, the time required for UAV to completely cover the desired area, increases. Here, the 
time that the UAV needs to provide the full coverage for the area by visiting all the stop points, 
is called delay. Hence, the delay depends on the travel time of the UAV between the stop points, 
and the time that UAV spends at each stop point for transmissions. Thus, the delay can be written 
as: 

r = Ttr + MTs ( 28 ) 


where Ttr is the total UAV travel time, M is the number of stop points, and Tg is the time that 
the UAV stays at each stop point. Clearly, the travel time depends on the travel distance and 
location of the stop points, and the speed of the UAV. The total travel time will clearly increase 
as the number of stop points increases. However, in general, the exact relationship between Ttr 
and M strongly depends on the locations of the stop points which do not necessary follow a 
fixed path/distribution for different values of M. As an example, it can be shown that the exact 
travel time for M = 3 and M = 4 is and ^ respectively, where v is the speed of the UAV, 
and Rc is the radius of the desired area. The residence time, Tg, depends on the multiple access 
method. If the UAV adopts a time division multiple access (TDMA) technique, the residence 
time will be a function of the number of stop points. In fact, a higher number of stop points 
corresponds to a smaller coverage region of the UAV. In this case, at each residence point, the 
UAV needs to provide service for a fewer number of users. Therefore, by increasing the number 
of stop points, the residence time can be decreased in the TDMA case. Considering a uniform 
distribution of the users, the residence time is approximately computed as: 




(29) 


where i is the service time of UAV for each downlink user, U is the number of downlink 
users, i?min is the coverage radius of the UAV which depends on M, the number of the stop 
points, and Rc is the radius of the desired area. However, if the UAV uses a frequency division 
multiple access (FDMA) technique, the users can be served simultaneously. In other words, the 
UAV does not need to use different time slots to serve the users. Therefore, if users are of 
homogeneous traffic type, the residence time of the UAV at each stop point does not depend 
on the number of the users, and hence it can be fixed. In this case, the residence time at each 
stop point will be constant and it does not depend on the coverage radius of the UAV and the 
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Figure 3: Total delay inereases as the number of stop points. 


number of stop points. As a result, i. In our model, we have eonsidered FDMA for 

multiple access. Hence, the residence time is the same for all values of M. In Figure we 
have shown the total delay versus the number of stop points for two values of residence time, 
and V = 10 m/s. As expected, the total delay increases as the number of stop points increases. 
Moreover, when the residence time of the UAV at each stop point increases, the additional delay 
due to a higher number of stop points increases. As we can see from Figure]^ for i = 20 s, the 
delay increases from 230 s to 480 s if the number of stop points increase from 3 to 10. However, 
for Ts 1 = 40 s the delay increases from 295 s to 690 s. Clearly, the delay and the number of stop 
points are directly related. It should be noted that, for our simulations, we consider the number 
of stop points as delay. 

Next, we derive the overall outage probability for a typical D2D user in the M time instances 
for the mobile UAV case. The outage probability is the probability of having at least one failure 
during M retransmissions. Assume that the relative location of the stop point with respect 
to the D2D user is {ri,hi) where r* is the distance between the projection of the UAV on the 
ground and D2D user and hi is the UAV altitude. Clearly, the distance between the user and 
UAV is \Xu^i\ = \/hj + ri^. For different time slots, the Rayleigh fading changes and can be 


considered independent [231. However, since the location of the D2D users do not significantly 
change during the multiple time slots, the interference from the D2D users are correlated. Then, 
the overall outage probability for D2D users can be found in the next theorem. 


Theorem 3. The overall outage probability for D2D users in M retransmissions considering 
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moving UAV is given by: 


Pout,d — 1 6Xp I 


IR'^ 


M 


1 - 


1 


1 + 


0\x\ 


M 


dec 


i=l 


exp 


Pd 


exp 




I 5 


(30) 


where M is the number of retransmissions, Iu,i is the interference from the UAV at retrans¬ 
mission, and Ei^f.) is: 


E; 


exp 


= Pt 


LoS,i exp ( -pr- I + -TNLoS,! 


jexp 


Pd 


Proof See Appendix C. ■ 

From Theorem 3, we ean observe that, inereasing M leads to a higher outage probability. 
In faet, as the number of stop points inereases, the UAV ereates a stronger interferenee on the 
D2D users. Consequently, Pout,d tends to 1 for M —>■ cx). However, the higher number of stop 
points for UAV enhanees the eoverage performanee of the downlink users. Henee, a tradeoff 
between eoverage performanee of downlink users and the outage of D2D eommunieations should 
be taking into aeeount. Moreover, Theorem 3 shows that, in order to guarantee that the outage 
probability does not exeeed a speeified threshold for different values of M, we should adaptively 
reduee the distanee between the D2D transmitter and reeeiver (do), or have orthogonal speetrum. 


V. Simulation Results and Analysis 
A. The static UAV scenario 

First, we eompare our analytieal results of the eoverage probabilities using numerieal simula¬ 
tions. Table II lists parameters used in the simulation and statistieal analysis. These parameters 
are set based on typieal values sueh as in Q and p0[ . Here, we will analyze the impaet of 
the various parameters sueh as the UAV altitude, D2D density, and SINK threshold on the 
performanee evaluation metries. 

In Figures and we show, respeetively, the D2D eoverage probability and the lower and 
upper bounds for the DU eoverage probability for different SINK deteetion threshold values. 
From these figures, we ean elearly see that, the analytieal and simulation results for D2D mateh 
perfeetly and the analytieal bounds for DU eoverage probability and the exaet simulation results 
are elose. Figures and show that, by inereasing the threshold, the eoverage probability for 
D2D users and DUs will deerease. 
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Table II: Simulation parameters. 


Description 

Parameter 

Value 

UAV transmit power 

Pu 

5 W 

D2D transmit power 

Pd 

100 mW 

Path loss coefficient 

K 

-30 dB 

Path loss exponent for UAV-user link 

O-d 

2 

Path loss exponent for D2D link 

CXu 

3 

Noise power 

N 

-120 dBm 

Bandwidth 

W 

1 MHz 

D2D pair fixed distance 

do 

20 m 

Excessive attenuation factor for NLoS 

T/ 

20 dB 

Parameters for dense urban environment 

B, C 

0.136, 11.95 



Figure 4: D2D coverage probability vs. SINK threshold 



Figure illustrates the system sum-rate (Gbps) versus the threshold for 1 MHz transmission 
bandwidth, X^u = 10“^, h = 500 m, and two different values of A^. By inspecting (23) in 
Section III, we can see that the rate depends on the coverage probability, which is a decreasing 
function of the threshold, /3, and an increasing logarithmic function of it. Clearly, for high values 













































20 



of (3, the received SINK cannot exceed the threshold and, thus, the coverage probabilities tend 
to zero. On the other hand, according to (20) and (21), as (3 increases, log 2 (l + (3) increases 
accordingly. However, since the coverage probability exponentially decreases but log 2 (l + (3) 
increases logarithmically, the average rate tends to zero for the high values of (3. Furthermore, 
for /3 —)■ 0, since log 2 (l + 13) tends to zero and the coverage probabilities approach one, the rate 
becomes zero. 

Figure |7] shows the impact of D2D density on the sum-rate. In this figure, we can see that 
a low D2D density yields low interference. However, naturally, decreasing the number of D2D 
users in an area will also decrease the sum-rate. For high D2D density, high interference reduces 
the coverage probability and consequently the data rate for each user. However, since the sum- 
rate is directly proportional to the number of D2D users, increasing the D2D density can also 
improve the sum-rate. According to the Figure 6, as the density of downlink users increases, the 
optimal Arf that maximizes the sum-rate decreases. This is due to the fact that, as \du increases, 
the contribution of DUs in the system sum-rate increases and hence increasing the rate of each 
DU enhances the system sum-rate. To increase the rate of a DU, the number of D2D users as 
the interference source for DUs should be reduced. As a result, the optimal decreases as as 
\du increases. For instance as shown in the figure, by increasing Xdu from 10“^ to 4 x 10“"^, the 
optimal Xd decreases from 0.9 x 10“^ to 0.3 x 10“^. 

It is important to note that the value of the fixed distance, d^, between the D2D pair signif¬ 
icantly impacts the rate performance. Figure shows the Usum as a function of the density of 
D2D users and do- From this figure, we can see that, the rate increases as the fixed distance 
between a D2D receiver and its corresponding transmitter decreases. Moreover, the optimal D2D 
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D2D density {D2D/m^) xio* 

Figure 7: System sum-rate vs. D2D density (number of D2D pairs per m^). 





Figure 8: System sum-rate vs. D2D density and do- 


density whieh leads to a maximum Csum, inereases by deereasing do. In faet, for lower values 
of D we ean have more D2D users in the network. For instanee, by redueing do from 8 m to 5 
m, the optimum average number of D2D users inereases by a faetor of 3. 

Figure shows the eoverage probability for DUs and D2D users as a funetion of the UAV 
altitude. From the DUs’ perspeetive, the UAV should be at an optimal altitude sueh that it ean 
provide a maximum eoverage. In faet, the UAV should not position itself at very low altitudes, 
due to high shadowing and a low probability of LoS eonneetions towards the DUs. On the 
other hand, at very high altitudes, LoS links exist with a high probability but the large distanee 
between UAV and DUs results in a high the path loss. As shown in Figure]^ for h = 500 m the 
DU eoverage probability is maximized. Note that from a D2D user perspeetive, the UAV ereates 
interferenee on the D2D reeeiver. Therefore, D2D users prefer the UAV to be at an altitude for 
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Figure 9: Coverage probability vs. UAV altitude. 


whieh it provides a minimum eoverage radius. As seen in Figure 8, for h ^ oo, the D2D users 
achieve the maximum performance. However, h = 800 m results in a minimum D2D coverage 
probability due the high interference from the UAV. 


Figure 10 shows the optimal UAV altitude that maximizes DU coverage probability versus the 
D2D users’ density. As we can see from Figure the optimal UAV altitude for downlink users 
decreases as the number of D2D users increases. This is due to the fact that a higher density of 
D2D users creates higher interference on the downlink users, and consequently the UAV reduces 
its altitude to improve SINK value for the downlink users. In other words, the UAV positions 
itself closer to the downlink users to cope with the high interference caused by the increased 
number of D2D users. From Figure we can see that, the optimal UAV altitude is independent 
of the fixed distance, do, between the D2D transmitter and receiver pair. In fact, the distance 
between D2D users does not affect the amount interference generated on the downlink users. 
Therefore, the optimal altitude of the UAV does not change if do changes. 


Figure 11 shows Usum versus the UAV altitude for different values of the fixed distance, do, 
the fixed distance between a D2D transmitter/receiver pair. The optimum values for the height 
which lead to a maximum Usum are around 300 m, 350 m, and 400 m for do = 20 m, 25 m 
and 30 m. Note that the optimal h that maximizes the sum-rate depends on the density of 
DU and D2D users. From Figure considering do = 20 m as an example, we can see that 

for h > 1300 m, the system sum-rate starts increasing. This stems from the fact that the DU 
coverage probability tends to zero and, thus, only D2D users impact Usum- Hence, as the UAV 
moves up in altitude, the interference on D2D users decreases and Ud increases. Moreover, for 
300 m < h < 1300 m. Figure shows that the coverage probability and, consequently, the 
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Figure 10: Optimal UAV altitude vs. D2D density. 



Figure 11: System sum-rate vs. UAV altitude. 


average rate for the downlink users decrease as the altitude increases. However, increasing the 
UAV altitude reduces the interference on the D2D users and improves the average rate for D2D 
users. In addition, in this range of h, since DUs have more contributions on Usum than the D2D 
users, Csum is a decreasing function of altitude. 


B. The mobile UAV scenario 

Next, we study the mobile UAV scenario. In this case, we can satisfy the coverage requirement 
for all the DUs. In fact, the UAV moves over the target area and attempts to serve the DUs at 
the stop points to guarantee that all the DUs will be in its coverage radius. 


Figure 12 shows the coverage radius of the mobile UAV when it is located at the optimal 
altitude as the D2D density varies. As expected, the coverage radius decreases as the D2D density 
increases. For instance, for £ = 0.6, when increases from 10“® to 10“^, the coverage radius 
decreases from 1600 m to 300 m. Moreover, by reducing the minimum coverage requirement of 
DUs, the UAV can cover a larger area. For instance, reducing e from 0.6 to 0.4 increases the 
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D2D density {D2D/m^) 


Figure 12: Maximum UAV coverage radius vs. D2D density (number of D2D pairs per m^). 


UAV coverage radius from 290 m to 380 m for = 10“^. Note that, since the main goal of 
the UAV is to provide coverage for the entire target area, to compensate for the low coverage 
radius, we should increase the number of stop points for serving the DUs and consequently a 
longer time is required for the full coverage. 


In Figure 13 , we show the minimum number of stop points as a function of the D2D user 
density. In this figure, we can see that, as expected, the number of stop points must increase 
when the density of D2D users increases. In fact, to overcome the higher interference caused by 
increasing the number of D2D users, the UAV will need more stop points to satisfy the DUs’ 
coverage constraints. For instance, when Xd increases from 0.2 x 10-^ to 0.8 X 10 the number 
of stop points must be increased from 3 to 8. Note that, when computing the minimum number 
of stop points for each Xd, we considered optimal values for the UAV altitude such that it can 
provide a maximum coverage for the DUs. Therefore, the UAVs altitude changes according to 
the D2D density. Moreover, as seen from Figure the minimum number of stop points remains 
constant for a range of Xd- This is due to the fact that the number of stop points is an integer 
and hence, for different values of Xd, the integer value will be the same. However, although the 
minimum number of stop points for two different D2D densities are the same, the UAV can 
transmit with lower power in the case of lower D2D density. 

In Figure we show the minimum number of stop points as a function of the UAV altitude 
for Xd = 10“^. Figure [14] shows that, for some values of h which correspond to the optimal UAV 
altitude, the minimum number of stop points is minimized. For example, the range of optimal 
/i for e = 0.4 and e = 0.6 is, respectively, 400 m < h < 500 m and 300 m < h < 350 m. As 
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Figure 13: Number of stop points vs. D2D density. 



Figure 14: Minimum number of stop points vs. UAV altitude. 

expeeted, the minimum number of stop points is lower for the lower value of e. 

Figure shows the tradeoff between the downlink eoverage probability and the delay whieh 
is eonsidered to be proportional to the number of stop points. In Figure we ean see that, in 
order to guarantee a higher eoverage probability for DUs, the UAV should stop at more loeations. 
As observed in this Figure, for = lO”"^, to inerease the DU eoverage probability from 0.4 
to 0.7, the number of stop points should inerease from 5 to 23. For a higher number of stop 
points, the UAV is eloser to the DUs and, thus, it has a higher ehanee of LoS. However, on 
the average, a DU should wait for a longer time to be eovered by the UAV that reaehes its 
vieinity. In addition, as the density of D2D users inereases, the number of stop points (delay) 
inereases espeeially when a high eoverage probability for DUs must be satisfied. For instanee, 
if Xd inereases from 0.5 x 10“^ to 10“"^, or equivalently from 50 to 100 for the given area, the 
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Figure 15: Minimum number of stop points vs. eoverage probability (eoverage-delay tradeoff) 



Figure 16: Overall D2D outage probability vs. number of retransmissions. 

number of stop points should increase from 4 to 9 to satisfy a 0.5 DU coverage probability, and 
from 20 to 55 for a 0.8 coverage requirement. 

Figure shows the overall outage probability for D2D users versus the number of retrans¬ 
missions. As the number of retransmissions (time slots) increases, the overall outage probability 
also increases. In other words, for higher number of time slots, the possibility that a failure 
happens during retransmissions, increases. Furthermore, since the UAV is an interference source 
for the D2D users, the higher number of stop points leads to a higher outage probability. From 
Figure we can see that, the increase in the outage probability of D2D users due to the UAV 
is 0.20 for M = 3, and is 0.38 for M = 7. Therefore, when the number of stop points increases 
due to the higher density of D2D users or a higher coverage requirement of the downlink users, 
the D2D communications are more prone to a failure. 
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VI. Conclusions 

In this paper, we have studied the performance of a UAV that acts as a flying base station in 
an area in which users are engaged in D2D communication. We have considered two types of 
users: in the network: the downlink users served by the UAV and D2D users that communicate 
directly with one another. For both types, we have derived tractable expressions for the coverage 
probabilities as the main performance evaluation metrics. The results have shown that a maximum 
system sum-rate can be achieved if the UAV altitude is appropriately adjusted based on the D2D 
users density. In the mobile UAV scenario, using the disk covering problem, the entire target 
area (cell) can be completely covered by the UAV in a shortest time with a minimum required 
transmit power. Moreover, in this case, we have derived the overall outage probability for D2D 
users, and showed that the outage probability increases as the number of stop point increases. 
Finally, we have analyzed the tradeoff between the coverage and the time required for covering 
the entire target area (delay) by the mobile UAV. The results show that, the number of stop 
points must be significantly increased as the minimum coverage requirement for DUs increases. 
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A. Proof of Theorem 1 
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where g is an exponential random variable with a mean value of one (i.e. g ~ exp(l)), (a) 
follows from the exponential distribution of g based on the Rayleigh fading assumption, and 
taking the expectation over and I2 (as random variables). Step (b) comes from the fact that 

and 12 are independent because the interference stems from different sources which are spatially 
uncorrelated. 

Here, E/^ and E/c are given by: 
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where the details of step (a) follow direetly from the results in [231. 
Finally, using (31), (32) and (33) Theorem 1 is proved. 
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B. Proof of Theorem 2 

The eoverage probability for a eellular user loeated at (r, (p) is written as: 
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Note that, there is no elosed-form expression for the eumulative distribution funetion (CDF) 


of the interferenee from D2D users [29| and Here, we provide lower and upper bounds for 
the CDF of interferenee. First, we divide the interfering D2D transmitters into two subsets: 

^l = {^B\Pdd^-^^g^>T}, 

$2 = {^^\Pdd-^^g, < T}, 

where T is a threshold whieh is used to derive the CDF of the interferenee from D2D users. 
Now, eonsidering the interferenee power from D2D users loeated in <hi and $2 as and 
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, we have: 
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(36) 


where (a) and (b) eome from the Rayleigh fading assumption and PGFL of the PPR 
The upper bound is derived as follows: 

F[Id<T] = l-F[Id> T] 

= l-{F[Id> T\Id,^, > T] P [Id,^, >T] + F[Id> T\Id,^, < T] P [Id,<s>, < T]) 
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= 1 - (P > T] + P [J, > < T] P < T]) 

= 1 - (l - P [$1 = 0] + P [Jrf > T\h,^, < T] P [$i = 0] 

= P[$i = 0] fl-P[/<i>T|$i = 0] 


( 37 ) 


Also, 
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(38) 


where (a) is based on the Markov’s inequality which is stated as follows: for any non-negative 
integrable random variable X and positive T, P{X > L) < Also, 1(.) is the indicator 
function which can only be equal to 1 or 0. Hence, the lower (T/) and upper (Uj) bounds for 
the CDF of interference become: 


Li{T)= 1- 


27rA,r(l + 2/a,)/T 


Oid 2 


\pd 
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exp I -7rAd( — ) r(l + 2/ad) ) , (39) 
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Ui{T) = exp ( -7rAd( ^ 


— ‘ijo-d 


r(l + 2 / ad) I . 


(40) 


Thus, we have Lj(T) < P{Id <T}< Ui(T). 

Finally, considering (34), (39), and (40), the lower bound and upper bound of the average 
coverage probability for DUs in the cell is expressed as: 
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+ / Wr)L, - iv) |dr. (41) 

PL.M = / PLos{r)U, - /v) ^dr 

+ J P^L,s{r)U, _ A,^ ^dr, ( 42 ) 

and Theorem 2 is proved. 


C. Proof of Theorem 3 

Consider 7 ^ j and gt, respeetively, the SINK and the ehannel gain (with exponential distribution) 
at retransmission, for 1 < i < M. The outage probability is the probability of having at least 
one failure during M retransmissions. Then, we have: 

Pout4 =1 - P [7d,i > I3, > f3] 
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(43) 


where (a) follows the assumption that the fading is independent in different retransmissions, and 
step ( 6 ) eomes from the faet that interferenee due to D2D users, interferenee from UAV, and 
noise are all independent. Also, 
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(44) 


where details of (c) ean be found in [23| where the eorrelation between D2D interferenee in 
different retransmissions is taken into aeeount. Finally, 
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-Pd'^^rtPAXuA-"^ 


Pr, 


(45) 


where step {d) is based on the fact that the interference from the UAV can be treated as 
independent in different retransmissions. 

Finally, using (40), (41), and (42), Theorem 3 is proved. 
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